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Magnetic anisotropies in quantum dots (QDs) doped with magnetic ions are discussed in terms of 
two frameworks: anisotropic g-factors and magnetocrystalline anisotropy energy. It is shown that 
even a simple model of zinc-blende p-doped QDs displays a rich diagram of magnetic anisotropies 
in the QD parameter space. Tuning the confinement allows to control magnetic easy axes in QDs 
in ways not available for the better-studied bulk. 
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I. INTRODUCTION 

Once the origin of magnetic ordering in a specific ma- 
terial is understood, it is often important to determine 
its magnetic anisotropy (MA) and hard and easy mag- 
netic axes in particular. A shift of focus towards MA has 
already occurred for the studies of bulk dilute magnetic 
semiconductors (DMS)^i^ but not yet fully for magnetic 
quantum dots (QDs) where it could play certain role, for 
example, in context of transport phenomena)^ the forma- 
tion of robust magnetic polaronsr'"— control of magnetic 
ordering)^"— nonvolatile memory;^ and quantum bitsJ^ 

In epilayers of (Ga,Mn)As, a prototypical DMS, the 
magnetocrystalline anisotropy energy (MAE) has been 
found to be a significant and often dominant source of 
MAii^"— caused by a strong spin-orbit (SO) coupling. It 
turns out that the easy axis direction depends on hole 
concentration, magnetic doping level as well as on other 
parameters. For example, when (Ga,Mn)As was used 
as a spin injector, the effects of strain (by altering the 
choice of a substrate) were responsible for changing the 
in-plane to out-plane easy axisJ^ While the strong SO 
couplingi^ is also present in p-type QD of zinc-blende ma- 
terials doped with Mn, its effect on magnetic anisotropies 
will be significantly modified by the confinement. The en- 
ergy levels in such 'nanomagnets,— where the Mn-Mn 
interaction is mediated by carriers, depend on the magne- 
tization direction bm = {nx,ny,nz)- It is often assumed 
that the interaction of magnetic moments with holes in 
quantum wells (QWs) or, equivalently in flat QDs, is ef- 
fectively Ising-like Here we quantify this assumption 
and explore MA using two frameworks: (i) an effective 
two-level Hamiltonian with a carrier g-tensor;^ which 
is widely employed also in theory of electron spin reso- 
nance, and (ii) MAE, which is commonly used to study 
bulk magnets. 

While previous studies focused on specific nonmagnetic 
QDsS^ and properties sensitive to system details (such as 
precise position of magnetic ions^^i^l), we explore more 
generic magnetic QD models, which can also serve as 
a starting point for more elaborate work. We consider 
a Hamiltonian comprising non-magnetic and magnetic 



parts, 

H^Hqd+Hcx- (1) 

The former encodes both QD confinement and SO inter- 
action, which is prerequisite for magnetic anisotropies, 
the latter expresses the kinetic-exchange coupling be- 
tween holes and localized magnetic moments. For 
transparency, we disregard the magnetostatic shape 
anisotropyS^ and assume that the QD contains a fixed 
number of carriers. We mostly focus on the case of 
a single hole; realistically, such system can be a II- VI 
colloidal^ or epitaxial. QD with a photoinduced carrier. 
Magnetic moments of the Mn atoms are taken to be per- 
fectly ordered (collinear) and are treated at a mean-field 
level. The magnetic easy axis is then the direction em for 
which the zero-temperature free energy F{eM) is mini- 
mized. In this article, we take two different points of view 
on F(eAf). On one hand, we discuss the lowest terms of 
F^bm) expanded in powers of the direction cosines of 
magnetization (n^ + riy + n1 = 1), inspired by the stan- 
dard 'bulk MAE phenomenology' and pay special atten- 
tion to the case of perfectly cubic QDs, F{eM) ~ Fo(ej\/). 
The anisotropies in Fq stem purely from the crystalline 
zinc-blende lattice. On the other hand, F{ei\i) acquires 
additional terms in systems with less symmetric confine- 
ment. We therefore discuss the anisotropic factors as a 
useful framework to handle such systems, e.g. cuboid 
QDs (orthogonal parallelepiped; extremal cases are a 
cube and an infinitely thin slab, i.e., a QW) and show 
how the expansion 

F{eM) ^ FoieAi) + AFiieAi) + A^F2{eM) + ■ ■ . (2) 

can be constructed using powers of A which reflects the 
anisotropy in ^-factors. Wc begin by discussing this lat- 
ter topic in Section II (quantity A is defined by Eq. ([5]) 
at the end of Sec. HA), then proceed to the phenomeno- 
logic (symmetry-based) expansions of Fq in Section III 
and conclude that Section with calculations of Fi in sit- 
uations that are beyond the applicability of the g-factor 
framework. 
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II. EFFECTIVE TWO-LEVEL HAMILTONIAN 

Since Hqd is invariant upon time reversal, its spec- 
trum consists of Kramers doublctsj^ To study the 
ground-state energy in the presence of magnetic mo- 
ments, we examine how these doublets are split by 
Hexi^Ai) where ej\/ is treated as an external parame- 
ter (related to classical magnetization; single-Mn doped 
QDs where the Mn magnetic moment behaves quantum- 
mechanically^ require different treatment) and represent 
them by an effective two-level Hamiltonian of Eq. ^ . We 
consider two example systems: a simple four-level one 
where completely analytical treatment is possible, and a 
more realistic envelope-function based model of a cuboid 
QD. 



A. Four level model 

Related to the Kohn-Luttinger Hamiltonian of a 
QWi^i^ the arguably simplest non-trivial model de- 
scribing anisotropy of a flat QD is 



.„ 1 
Hi = aJ^ + -heM ■ J 



(3) 



representing hole levels in a zinc-blendc structure 
whose confinement anisotropy and exchange splitting are 
parametrized by a and /i, respectively (the term aJ^ 
implies that the strongest confinement is along the z- 
direction and this term also encodes information about 
the SO coupling). Jx,y,z are 4x4 spin-§ matrices. In 
terms of Eq. ([T]), we now choose H ~ Hi and the first 
(second) term in Eq. (jS]) plays the role of Hqd {Hex)- 
Anisotropic behavior of eigenvalues of iJi, to linear or- 
der in /i/a, is illustrated in Fig.[TJa). It can be extracted 
from the exact eigenvalues. 



— a ± —n 
4 6 

— a ± —n 
4 6 



-h^ =F -ah 
9 3 

-h^ =F -ah 
9 3 



(4) 
(5) 



in the C3jSG Tlx = 1 (or ny = I), shown in Fig.jljb), which 
clearly differ from the case = I where the eigenvalues 

9a/4±/i/2and 



are strictly linear functions of h {E^f^ 



a/4 ± h/6); subscripts refer to the i?^,j(0) = 9a/4 



^ih 

('heavy-hole', HH) and E^f^{0) = a/4 ('light-hole', LH) 
doublets, respectively. In the limit of weak exchange, 
h/a ^ 1, splitting of each of the Kramers doublets is 
symmetric and it can be characterized by three param- 
eters \dE/d{hnp)\, p = x,y,z, for /i — > as depicted 
in Fig. [Ija). These parameters can be plausibly called, 
by analogy with the Zecman effect, the anisotropic g- 
factors gp. From Eqs. ([I]),©, we straightforwardly ob- 
tain igx,9y,gz) = (0,0,1/2) and (1/3,1/3,1/6) for the 
HH and LH doublet of the Hamiltonian Hi , respectively. 
This result is known from the context of QWsi^i^ We 
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FIG. 1. (Color online) Splitting of levels E{h) in a flat QD 
described by Eq. (a) For the particular Kramers doublet, 
E{h) depends on eM and the factors (by convention non- 
negative) are dE/dh — {gx, gy, Qz)- (b) Beyond the linear 
regime in h/a, dE/d{hnx) will be different for the upper and 
lower level of the split doublet, it will depend on h and may 
even change sign, indicating that the Haft of Eq. ((6| based on 



emphasize that these ^-factors of the model specified by 
Eq. ^ are independent of the parameters a, h (except 
for the requirement h <^ a which represents the /i — > 
limit). 

If we focus on one particular Kramers doublet, it is 
straigtforward to show that Hi projects to 



H 



off 



[nxgxTx + n-ygyTy + n^g^T^l 



(6) 



for a suitably chosen basis \Ki), \K2) of the doublet. Here 
Ti are Pauli matrices and we have mapped two eigenstates 
of the original Hamiltonian Hqd on a pseudospin \ f\ = 
1/2 doublet \+), I-), where t,\±) = ±|±). For H = Hi, 
the eigenstates are only four-dimensional (spanned by 
the \.h = 3/2), |J, = -1/2), |J, = 1/2),|J, = -3/2) 
basis). We present another example of H in Sec. IIB 
where advantage of the projection becomes more appar- 
ent. The choice of basis |-|-), |— ) is crucial to obtain 
Hcff in the simple form ([5]); considering the HH dou- 
blet: 1+) = \J, = 3/2), I-) = IJ;, = -3/2) leads 
to Eq. ^ while for other basis choices the mapping 
Hex — {h/3)eM ■ J I— ^^off = heM ■ g ■ t may lead^ to 
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non-symmetric tensor g = gij, i,j £ {x, y, z}. In general, 
if the mapping is to produce gij = diag {gx,gy,gz) the 
'suitable choice of the basis \K-i), |if2)' where \Kx) i-^ |+) 
is such that {Ki\Jx,y\Ki) = 0, {Ki\J^\Ki) > (and l/fj) 
is the time- reversed image of \Ki) which is mapped to 

Let us now consider a general system described by 
Eq. H]) . Assuming that the downfolding of H into HeS is 
possible for given \Ki), \K2) (this assumption is discussed 
in Appendix A), the anisotropic factors can readily be 
determined as dE/dh for the particular Kramers doublet 
level E. This is equivalent to perturbatively evaluating 
the effect of Hex on two degenerate levels to the first 
order of h as follows: (i) specify the Kramers doublet 
of interest, and find any basis |A"i), \K2) of this dou- 
blet, (ii) extract the operators tx,y,z from Hex by taking 
tp = dHex/d{nph) (for example, tx = J^/S for HcS ap- 
pearing in Hi), (iii) evaluate their matrices 



{Kl\ix,y,z\Kl) {Kl\ix,yMK2) 
{K2\L,y,z\Ki) {K2\ix,y,AK2) 



(7) 



in the two-dimensional space spanned by \Ki), \K2), 
and (iv) the non-negative eigenvalue of tp equals gp 
(p = x,y,z). We emphasize that while gp depends on 
system parameters in Hqjj and Hex, it also depends on 
which Kramers doublet we choose. Higher doublets be- 
come relevant for QDs containing higher (odd) number 
of holes, for example. 

The effective Hamiltonian in Eq. (jB]) can be used for 
various purposes, e.g., for studies of fluctuations of mag- 
netization in magnetic QDs^, spin-selective tunneling 
through non-magnetic QDs^i or excitons in single-Mn 
doped QDsj^ If the magnetic easy axis is of interest, 
the (/-factors immediately provide the answer: F{eM) 
based on Eq. ([S]) is minimized for e^j in the direction 
of the largest gp (e.g. for the HH doublet in Fig. [IJa), 
it is = 1 because gz > gx,gy)- If the full form of 
F(ej\/) is needed (e.g, for ferromagnetic resonance^), it 
can be straightforwardly obtained by diagonalizing the 
2x2 matrix of Hcs- Assuming gx = dy, the (modulus of 
the) eigenvalue can be expanded in terms of parameters 
A and k as derived in Appendix B. It is meaningful to 
call 



A={9l-9l)/{9l+9l) 



(8) 



the asymmetry parameter since it vanishes in a perfectly 
cubic QD {gx = gy = gz) and it is with respect to this 
parameter that we can identify 



AFi{BM) = -Aknl 

A^F2{eM) = +\A^k{2nl - 1)^ 
o 

in Eq. ^ to linear order oi k oc h. 



(9) 
(10) 



B. A cuboid quantum dot model 

With this general scheme at hand, we take one step 
in the hierarchy of models towards a more realistic de- 
scription of magnetic QDs. We consider a zinc-blendc 
structure p-doped semiconductor shaped into a cuboid 
of size Lx X Ly X Lz such as can be described by four- 
band Kohn-Luttinger Hamiltonian.— Also in this system, 
H = H2 is a sum of Hex and Hqd but this time, i?QD 
comprises of blocks {mxmymz\HKL\^'x'^'y^'z) with 



Hkl — 



{(71 + ii^v - '^12[jipI + Jivi + Jim 



2mQ 

-^-faliJxJy + JyJx)P xPy 



c.p. 



(11) 



Here, \mxmymz) denotes the basis of envelope functions, 
7i,2.3 the Luttinger parameters, mg the electron vacuum 
mass, Px.y.z the momentum operators and c.p. denotes 
the cyclic permutation (see Appendix C for details) . The 
envelope function is conveniently developed into har- 
monic functions with rup — 1 nodes in the p ~ x,y, z 
direction: 



{r rnxmyniz 
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iVsin^^sin^^sin^— . (12) 
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We have introduced the dimensionless aspect ratios 
Ax. J/ = Lx^y/L and the normalization factor N . Our sys- 
tem can be viewed as an infinitely deep potential well 
with V{x,y,z) = for < a; < L^, < y < Ly and 
< z < Lz = L and infinite otherwise. 

For fixed material parameters (Luttinger parameters 
in ratios 72/71, 73/72) and QD shape (A^, Aj^), all matrix 
elements of all blocks {mxmymz\Hji]^\m'^m'ym'^) scale 
as 1/L^. The spectrum, consisting of Kramers dou- 
blets which occasionally combine into larger multiplcts, 
is specified by a sequence of dimensionless numbers E/£q 
where 



£0 = h^n^7i/i2moL''). 



(13) 



For a cubic QD [A^ = Ay = 1; see Fig. [D(a)] the s-like 
state shown in the inset of Fig. [2][a) forms a quadruplet, 
and depending on the value of 72/71 (and to somehow 
lesser extent also of 73/72) this state competes with the 
next doublet for having the lowest energy. The critical 
value (see Appendix C) 



CR = (2 + 128/97r2 



0.29 



(14) 



can be taken to distinguish materials with small (72 /71 < 
cr, ground state quadruplet) and large (72/71 > c/j, 
ground state doublet) splitting between light and heavy 
holes in the bulk; these can be ZnSe and CdTe, respec- 
tively, their values of 72/71 based on approximating 72 
and 73 by their average 72 = (72 + 73)/2 are indicated 
in Fig. [Ha). By numerical diagonalization we have de- 
termined the lowest 7 Kramers doublets in slightly de- 
formed QDs {Xx = Xy = X = 1.01) in these materi- 
als (71/2/3 = 4.8/0.67/1.53 for ZnSc and 4.1/1.1/1.6 for 
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CdTc)2§ and executed the procedure (i)-(iv) above to ob- 
tain the g-factors which are hsted in the table on the right 
of Fig. [2] {gx — gy due to A^; = A^). To avoid confusion, 
we remark that in (i), \Ki), {^2} are vectors of dimension 
864 in the basis \mxmymz) (g) | J^) (see the discussion of 
cut-off in Appendix C) and in (ii), tx = {l/3)Jx (81 l^yz, 
where Ti-xyz is the identity operator in the space of the 
envelope functions given by Eq. p2p . Evaluation and di- 
agonalization of the 2x2 matrices in Eq. ([7]) requested in 
(iiijiv) is performed numerically. The possibility to map 
the action of Hex = (/i/3)ej\/ • J (g) Tlxyz on the Kramers 
doublets \Ki), |i^2) implied by Hqd of a cuboid p-doped 
QD is discussed in Appendix A. 

The slight deformation of the QD makes the quadru- 
plet split into two doublets (with energies 71.7 and 
71.9 meV for ZnSe) whose (7-factors approach (0,0,1/2) 
and (1/3, 1/3, 1/6). Similar situation occurs for the dou- 
blet pair with energies 52.8 and 53.0 meV for CdTe. The 
actual ground state in this material is, however, a dou- 
blet of different orbital character than the quadruplet 
(we stress that this is due to the confinement, see Ap- 
pendix C); it evolves from the E = 6£o level of 72/71 = 
as shown by the solid line in Fig. ^a) and its ^-factors 
are isotropic, (1/6,1/6,1/6) in the limit A — >■ 1. This 
doublet, however, remains the ground state only in rather 
symmetric QDs (A w 1.25 in CdTe) and for more strongly 
deformed QDs, the lower doublet of the E = 3£o (at 
72/71 = 0) quadruplet becomes the ground state just 
as it is the case for ZnSe for arbitrarily small deforma- 
tions A > 1. In Fig. mb), we show how the factors 
of the CdTe QD ground state depend on A beyond the 
mentioned value ~ 1.25. These results, including the g- 
factors, are independent of the QD size L, except for the 
energies which scale as 1/L^ as mentioned above. 

From Fig. one may conclude that the Ising-like 
Hamiltonian is often an excellent approximation {gx — 
gy = 0, as others assumei^i2ii2^"i^) for the lowest 
Kramers doublet. To be more specific, we now discuss 
materials with small and large HH/LH splitting sepa- 
rately. For 72/71 < Cfl, the out-of-plane (7-factor (gz) 
overwhelmingly exceeds the in-plane one {gx = gy) even 
for minute deformation of the QD; this can be seen from 
the numeric ZnSe data in Fig. [5] We find gz — 0.464 
and gx = 9y = 0.012 for A — 1 as small as 0.01. For 
CdTe, which represents the other class (72/71 > cr), we 
find similar values {gz = 0.418) for the second Kramers 
doublet while the lowest doublet remains rather isotropic 
{9x = 9y = 0.166 and g, = 0.164). As we make the QD 
deformation larger, these two doublets cross, so that the 
ground state doublet is Ising like while the second lowest 
doublet remains more isotropic. This crossing occurs for 
A « 1.25 in CdTe and data in Fig. ^h) are only shown 
for A > 1.25. 

We now elaborate on the properties of the low-energy 
sector of H2 (at h = 0). Coupling between blocks of dif- 
ferent Inixmyinz) vanishes when 73/71,72/71 — > 0, and 
Eq. ^ becomes in this limit the exact effective Hamilto- 
nian of the lowest four levels {mp — 1 for all p = x, y, 2). 
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FIG. 2. (Color online) (a) Levels in a cubic dot (with 73 — 72) 
in units of £0 defined by Eq. (|13l) . Solid lines indicate ana- 
lytic result obtained when mixing between remote levels is 
disregarded. Note that their crossing (which we use to dis- 
cern the weak and strong HH/LH splitting materials, dashed 
line) is very close to the actual crossing when level mixing is 
taken into account. Values representing ZnSe (72/71 ~ 0.23) 
and CdTe (72/71 ~ 0.33) QDs are indicated. Inset: squared 
wavefunction modulus of the ZnSe QD ground state in the 
z = L/2 section, (b) Dependence of the factors associated 
with the ground state Kramers doublet in a CdTe QD on 
its shape (A^ — Xy = X). Right: Energies and g-factors in 
slightly deformed QDs (A = 1.01) for the lowest 7 Kramers 
doublets for ZnSe and CdTe, where £q « 28 meV and 24 meV, 
respectively, for L = 8 nm. 



They form a quadruplet for A = 1, which splits into two 
doublets upon deformation of the QD; we can see it by 
writing 



(llllffj^Lllll) = 3£o 



l4/(A)-jf(l-A-2) 



272 
3 71 



(15) 

where I4 is a unit 4x4 matrix and /(A) is a certain 
function with lim^^i /(A) = 1. The lower doublet of this 
4x4 effective Hamiltonian has gz = 1/2 (when A > 1 and 
72 > 0) and therefore the values of gz deviating from 0.5 
(appearing in Fig. ^ occur only due to admixtures from 
higher-orbital {mp > 1) states of LH character. Indeed, 
going from ZnSe to CdTe, the mixing becomes stronger 
and gz of the HH-like level drops from 0.464 to 0.418 
(A = 1.01, numerical data in Fig. ^. While Eq. © 
may remain the effective Hamiltonian of the two doublets 
originating from \mxmymz) ~ |111) even for 72/71 > 
(CdTe levels of 52.8 and 53.0 meV in Fig.©, for A close 
to 1, there is the more isotropic doublet on the stage 
(49.9 meV in Fig. [2]). Nevertheless, if A is sufficiently 
large, the J| term in Eq. (|lip will eventually dominate, 
it will suppress all mixing between HH and LH states 
and the lowest doublet will again approach {gx, dy, dz) = 
(0, 0, 0.5) as it is shown in Fig. [2;b). 



5 



III. MAGNETOCRYSTALLINE ANISOTROPY 
ENERGY 

In analogy to the bulk systems, even cubic QDs retain 
anisotropies. However, these cannot be described within 
the previous framework: for instance, gx,gy,gz arc all 
equal to 1/6 in the cubic CdTc QD ground state hence 
A = in Eq. ([5]). One could replace gij by a higher 
rank tensor to capture these effects, but MAE formalism 
of bulk magnets seems more customary and informative. 
Unlike the (^-factors, MAE analysis does not invoke the 
concept of Kramers doublets. The zero-temperature free 
energy F(eM) of a magnetic QD with a single hole is 
now simply the lowest eigenvalue of Eq. ([1]) and it can be 
expanded in powers of rij. The lowest terms compatible 
with cubic symmetry are^S 

Fo = K,{ni +nl + n\) + 27K,2nlnlnl (16) 

For data calculated by numerically diagonalizing H = H2 
(model described in Sec. IIB) it turns out that Eq. 
suffices to obtain good fits; for instance, lower solid line 
in Fig. Ela) corresponds to Kc = 0.83 meV and Kc2 ~ 
0.075 meV with easy axis along [111]. There we have 
chosen Cdi_a;Mn3;Te as the material, L ~ \Q nm and 
/i = 50 meV which corresponds to h = JpdNMnSi\in 
with X « 2.3% (we take^ \ Jpd\ = 60 meV • nm^, Smu = 
5/2 and Nmu = ^x/af with CdTe lattice constant a; = 

0. 648 nm). Results in Fig.|3]are again subject to scaling, 
similar to the non-magnetic spectra in Fig. [21(a). When 
the material parameters (specifically, 72/71 and 73/72) 
are fixed, the spectrum of H2, expressed in the units of 
£0, depends on a single dimensionless parameter 

Z = h/£Q = 2?7io/iLV(7i7r2ft2). (17) 

This scahng relates the spectra of e.g. cubic dots of dif- 
ferent sizes and Mn contents (if their respective values 
of Z are equal). Data in Fig. [3] therefore apply both to 
X = 2.3% at L = 16 nm (if left as they are) and x = 9.2% 
at L = 8 nm (if scaled by a factor of 4). It turns out that 
the (7-factor analysis presented in the previous section is 
meaningful for Z < 0.1 while now wc have stepped out of 
this limit. When the exchange field h becomes stronger, 
levels cross and cease to depend linearly on h as required 
by Eq. (0; for F = ffi, this is illustrated in Fig. [^b). 
This limit was determined for CdTe cubic QDs but it will 
typically not be too different for other materials and/or 
aspect ratios A unless accidental (quasi) degeneracies oc- 
cur at Z = 0. 

MAE shown in Fig. [3] describe systems well beyond 
this limit of small Z (linear regime). Wc first focus on a 
perfectly cubic CdTe QD where there are no anisotropies 
in the linear regime. As already mentioned, the lowest 
energy hole state in Fig. Ela) exhibits a [111] easy axis 
with Kc = 0.83 meV at L = 16 nm and h = 50 meV, 

1. e. Z « 2.8; (this corresponds to a realistic x 2.3% 
Mn doping). In bulk DMSs, [111] would be an unusual 
magnetic easy axis dircctioi>i^ and we surmise that the 
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TABLE 1. Magnetic anisotropy constants (in meV) for a 16 x 
16 X 16 nm^ (cubic) and 16 x 16 x 8 nm^ (deformed) ZnSe 
and CdTe magnetic QD as a function of exchange splitting 
(or dimensionless parameter Z as for CdTe). 



reason for this is that for instance in (Ga,Mn)As grown 
on a GaAs substrate, there is a sizable compressive strain 
which prefers either parallel or perpendicular orientation 
of ejv/ with respect to the growth axis. 

Wc note that in a QD containing two holes (closed- 
shell systeroii) the anisotropies will also be present and 
they will be different from the single-hole case. Free en- 
ergy, taken as a sum, Fo(e7\/) = Ei + E2, of the lowest 
two single- hole states [shown e.g. in Fig. [31[a)], is not a 
constant independent of gm as one could naively expect. 
This intuition reflects H^s in Eq. © where the two hole 
states have opposite spin (hence their energies add up to 
zero). Once we leave the linear regime (Z > 0.1), i?eff 
ceases to be a good approximation. Qualitatively, the 
same behaviour is found for ZnSe (not shown) , a smaller 
value of Kc ~ 0.41 meV is accounted for by the smaller 
HH/LH splitting. The value of this constant is a com- 
plicated function of system parameters and it can even 
change sign as shown in Fig.[31[c) where Kc = —0.63 meV. 
Parameters used in this figure (71/72/73 = 4.0/1.5/1.6 
and h ~2Q meV) do not strictly correspond to published 
values of any semiconductor but they can be viewed as 
reasonable given the uncertainty in experimental deter- 
mination of the Luttinger parameters. Dependence of 
the anisotropy constants for ZnSe and CdTe QDs on h is 
summarized in Tab. HI 

Let us now return to non-cubic QDs. As already 
explained, the sizable (/-factor anisotropies shown in 
Fig. [HJb), relevant to the case of weak magnetism (Z ^ 
1), translate into an additional term AFi = Kunj, in the 
free energy of Eq. ^ where Ku = —kA up to linear or- 
der in Z (X k. Typically, K^ exceeds Kc already for small 
QD deformation (A slightly over one) and the data in 
Fig. EKb) imply Ku almost an order of magnitude larger 
than Kc for A = 2 (see also data in Fig. [2] where gz ^ gx)- 
Regardless of the contributions to Ku of higher order in 
Z, data in Tab. U imply an out-of-plane easy axis (in the 
[001] direction) as it is the case in QWs. However, upon 
deforming of a QD the easy axis does not abruptly jump 
from [111] to [001] but smoothly interpolates between 
these two directions. Similar effect, easy axis shifting as 
a function of some system parameter, is also known in 
bulk DMS [(Ga,Mn)As epilayers in particular, sec Fig. 8 
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FIG. 3. (Color online) Magnetocrystalline energy as a function of magnetization direction (Ei); the data labelled E2 are 
explained in the text. CdTe QD with 2.3% Mn (a) 16 x 16 x 16 nm^, (b) 16 x 16 x 8 nm^. (c) Fictitious material with 
parameters described in the text; note that the sign of Kc implied by Eq. (|16l) has changed compared to (a,b). (d) Color-coded 
easy axis positions for CdTe QDs as a function of aspect ratio (A) and effective exchange splitting Z. Black squares (1) indicate 
easy plane perpendicular to z-direction, hollow squares denote an isotropic magnet; white region (2) corresponds to easy-axis 
[001]; red squares (3) denote systems with [111] easy axis which gradually shifts towards [001] with increasing A. This plot is 
universal as far as L is concerned. 



in Ref. [l^ . Easy axes as a function of QD shape (oblate 
dots, A > 1) and effective exchange splitting Z are sum- 
marized in Fig. [3ljd) and the mentioned gradual shift of 
easy axis is indicated by shading between regions (3) and 
(2) (easy axes [111] and [001], respectively). On the other 
hand, the easy axis position changes abruptly between 
(1) and (3) or (1) and (2); region (1) corresponds to easy 
axis in the plane perpendicular to [001] (with anisotropies 
within this plane being very small) . The abrupt changes 
reflect ground state crossings, such as the one with A de- 
scribed below Eq. (|14p , while the gradual ones stem from 
level mixing caused by Hex- 

Finally, we comment on MA in QDs occupied by more 
than one hole. As already mentioned above, one pos- 
sible approach is to discuss open-shell and closed-shell 
systems separately. This notion is based on the concept 
of the QD being an artificial atom whose levels are or- 
ganized into shells comprising of spin-up and spin-down 
orbitals. Whenever a shell is completely filled (closed), 
the numbers of spin-up and spin-down carriers are equal 
hence their total spin is zero. If the QD is magnetically 
doped, no magnetic ordering is expected and also no MA. 
However, strong SO coupling puts this concept into ques- 
tion since it mixes different shells and also invalidates 
the spin-up and down labels of individual orbitals. The 
MA as a function of particle number Np strongly varies, 
both quantitatively and qualitatively. By comparing the 
Np = 1 and Np = 2 cases of a cubic CdTe QD, that is 



FoieAi) = El and Fo{eM) = Ei + E2 of Fig. ^a), we 
find that while the easy axis [111] in the former case is 
relatively 'soft' (energy difference between ejv/]][lll] and 
[110] is 'only' « 0.1 meV), the QD with two holes has a 
'robust' easy axis [110] and the corresponding minimum 
in Fo{eM) is as deep as 0.3 meV. MA as a function of 
Np displays rich behavior and one can therefore envision 
control of nanomagnetism by electrostatic gating, illumi- 
nation (used to photoinduce carriers) and possibly also 
temperature, known to alter the magnetic ordering in the 
bulk-like structures J^i^ 



IV. CONCLUSIONS 

We have discussed two approaches to magnetic 
anisotropies in quantum dots (QDs) described by a 
generic model in Eq. An effective Hamiltonian for 
individual Kramers doublets allows to express the ener- 
getics of a magnetically doped QD in terms of only three 
parameters (anisotropic g-factor) if the exchange split- 
ting due to the magnetic ions is relatively small. On 
the other hand, if the exchange splitting is large or the 
QD's symmetry is too high, the symmetry-based expan- 
sion of the magnetocrystalline energy in powers of the 
direction cosines of magnetization may in principle con- 
tain infinitely many terms (each of them quantified by 
one parameter). Focusing on manganese-doped semi- 
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conductor QDs, we find that only first few terms are 
appreciable, present their values and show in Fig. ^d) 
a diagram of magnetic anisotropics in the QD param- 
eter space. While we focus on a relatively small pa- 
rameter range in that diagram, and the barriers be- 
tween individual free energy minima are relatively low, 
it demonstrates that the QDs may have rich magnetic 
anisotropics. In spintronicS j^^'^^i^^ these systems could 
thus enable confinement-controlled multi-level logic. Our 
results provide a starting point for further studies of 
nanoscale magnetism in QDs. Such studies could relax 
the mean-field approximation, include multiple-carrier 
statesr^i^ or the effect of strain. 
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APPENDIX A 

The downfolding of Hqd + Hex to HcS is indeed possi- 
ble for the two example systems discussed in Sec. II A 
and IIB. To prove this, we first transform the basis 
\Ki),\K2) to \K[),\K!2) where of Eq. © is diagonal 
and then verify that the diagonal elements of and ty 
vanish. This procedure has to be applied to each Kramers 
doublet of interest. In the case of Hi in Eq. this is 
done simply by construction (e.g. \K[), j/Vj) for the up- 
per doublet in Fig. Ufa) is just \ J^ = 3/2), \J^ = -3/2)). 
In the model described by H2, one can split the Hilbert 
space into two disjunct subspaces Jifi, ,^1 and the above 
assertion can be shown to hold if |A'() G M\ and jATj) G 
3^1. (The decomposition M\ ® M\ relics on Hex being 
independent of space coordinates; relaxing the mean-field 
treatment of Mn magnetic moments thus introduces cor- 
rections to HcS-) Finally, one adjusts the relative phase 
between |Ar() and \K2), so that the matrix tx is real and 
ty purely imaginary. 



APPENDIX B 

This Appendix explains the relation between the 
anisotropic g- factors and Eq. The eigenvalues of Hcb 
are two numbers of equal magnitude and opposite sign, 
the lower of which is 

- h^nlg^+nlgl + nlgl (18) 

Let us consider for example single hole in a cuboid QD of 
dimensions XL x XL x L (such as it corresponds to data 
in Fig. so that gx = Qy Expression (jl8p which now 



equals F{eM) can be rewritten as 

\ gi + gi 

and developped in terms of a small parameter A = {g^ — 
9x)/(9z ~^ 9x) which quantifies the QD asymmetry as 

-fc(l - iA) - Aknl + \A^k{2nl - 1)^ + . . . (20) 

where k = hy^{g^+g^)/2. The first term does not de- 
pend on the magnetization direction, hence it can be dis- 
regarded for the purposes of magnetic anisotropy analy- 
sis. 

APPENDIX C 

We derive Eq. in this Appendix and discuss the 
details of the model considered in Sec. IIB. Energies E/Eo 
in Fig. [2][a) are calculated by numerical diagonalization 
of H2 with h ~ 0, a matrix constructed of 4 x 4 blocks 
{mxmymz\HKL\'rrL'^'m'yrn'^) / Eq introduced at the begin- 
ning of Sec. IIB. The basis of Hqd consists thus of direct 
product states \mxmymz) ® \ Jz) where \ Jz) are the four- 
spinors of total angular momentum J = 3/2 which are 
eigenstates to J^. For practical purposes, we cut-off the 
basis by m.x, my, ruz < 6, resulting in Hqd of dimension 
864. Eigenvalues are typically converged to better than 
0.1 mcV for this cut-off. 

The matrix Hqd/£q is block-diagonal for 72 = 73 = 
and the block mx,my,mz has a four- fold degenerate 
eigenvalue 

{mx/Xxf + {my/Xyf + ml. (21) 

Dimcnsionless energies on the left of Fig. [2][a) correspond 
to Aa; = Aj, = 1 and are hence integers. The lowest 
level E/Eq = 3 belongs to {mx,my,mz) = (1, 1, 1) while 
the first excited state E/£o = 6 entails an additional 
threefold geometric degeneracy corresponding to orbital 
states (1,1,2), (1,2,1) and (2,1,1); the E/Sq = 6 level 
for 72 = 73 = is thus twelve-fold degenerate. 

Next, wc can treat the HH-LH splitting as a perturba- 
tion when 72 and 73 arc turned on. In the lowest order, 
mixing between different {mx, my, mz) blocks can be ne- 
glected except for the case when their energies were equal 
at 72 = 73 = as in the case of the three blocks of the 
E/£o = 6 level. With coupling to the remote levels dis- 
regarded, we are left with a 12 x 12 matrix in this case 
which can be diagonalized analytically. It turns out to 
have two four-fold degenerate eigenvalues 

+ , 64 72 / n 8l7r4\ , ^ 

and two two-fold degenerate ones 

± 12872 / ^97r2\ 

^2/^0 = 6-— -(^.T^j. (23) 
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The lowest of these four energies is E2 and it is shown 
in Fig. [2Ia) for s = 73/72 = 1 as a soUd Une which 
crosses the horizontal line E/Eq = 3 corresponding to the 
(rnx,my,mz) = (1, 1, 1) quadruplet which does not shift 
in energy to the first order of this perturbation analysis. 



Eq. (|T4)) is the solution of E2 = Sfo for 72/71 under the 
assumption s = 1. Such level crossing (as a function of 
72/71) is genuinely due to the confinement and no level 
crossings occur in in bulk as long as < 72/71 < 1/2. 
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